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Abstract 

This text describes a generalization of the analytic signal taaboil.ll946 h ap- 
proach for the definition of instantaneous amplitude and phase to the case 
of multivariate signals. It was originally written as an appendix for an- 
^H i other paper, where the determination of the locally dominant oscillatory 

direction (the instantaneous amplitude) described here is used as a pre- 
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j^ ' processing step for another kind of data analysis. The text is reproduced 



cd 



C/3 






o 



^ 



in a 'standalone' form because the procedure might prove useful in other 



——k • contexts too, especially for the purpose of phase synchronization analysis 

(Rosenblum et al., 1996) between two (or more) multivariate sets of time 



series (Pascual-Marqui, 2007). 



C/3 . The local oscillatory behavior of a real-valued univariate signal x{t) is com- 

^~>' monly characterized using the corresponding complex-valued analytic signal 

z{t). It is obtained by combining x{t) with an imaginary part. 



z{t)=x{t) + iy{t), 



^ . which is defined as the Hilbert transform of x, 

en ■ 

where P. V. denotes the Cauchy principal value of the integral. Under the condi- 
tion that x{t) is dominated by a single frequency component, its instantaneous 
amplitude A{t) and phase ^{t) can be determined via the analytic signal ac- 
^*) i cording to 

A{t) = \z{t)\, ,/)(0 = argz(0, 

so that 

x{t) = A{t) cos(p{t) 



or 

z{t)=A{t) exp(i<^(0). 

The terms amplitude A and phase ^ as they are used here can be interpreted 
such that they specify the parameters of a strictly periodic sinusoidal oscilla- 
tion which locally matches the behavior of the observed signal x{t) at a given 
instant t. In particular, cj){t) attains the value (or equivalently, an integer mul- 
tiple of 2tz) whenever the actual value of x{t) coincides with the associated 
instantaneous amplitude A{t). 




Figure 1: Determination of local ellipse axes. The trajectory formed by the multivariate 
signal is locally matched to an elliptical orbit, which is defined by the data vector x at 
a given instant and the corresponding vector y from the signal's channel-wise Hilbert 
transform as conjugate axes. The main axes of the ellipse, a and h, are obtained using 
the associated multivariate instantaneous phase (p. 



These properties of the analytic signal can also be utilized to determine 
the parameters of the locally matching oscillation for a multivariate signal 
x{t) = (x,(t)) (f = 1 . . .K). We assume that each component signal x,(f) is 
dominated by a single frequency and that the frequencies of different signals 
are practically identical. Using y(f) to denote the channel-wise Hilbert trans- 
form of x{t) and z(f ) for its channel-wise completion to the analytic signal, the 
local extension of the signal's oscillatory behavior for instant t is obtained with 

zt{e)=z{t) exp(i0), 

parametrized by G [0,27r]. Its real part 

xt(e) =x{t)cose-y{t)sm.e 

gives the multivariate oscillation that locally matches the behavior of the signal 
at instant t; its trajectory is an elliptical orbit with conjugate axes defined by the 
vectors x{t) and y{t). 

From these conjugate axes, the main axes a{t) and b{t) of the local ellipse 
can be calculated (Fig. [1]|. It proves useful to do so via introducing a global 
(channel-independent) instantaneous phase (j>{t), such that for (p{t) G {0, ^tt, 
71, J n} (or equivalents), x{t) coincides with one of the main axis vectors or its 
negative. This is achieved choosing 

1 2x(t) -yit) 

(p{t) = -arctan- V ; i/v ; 



2 |f(f)|2-|^(f)|2- 

Since the resulting values in the range [—j,j] cover only one quarter of a 
cycle, the outcome may be transformed into an equivalent but more useful 



representation via a standard "unwrapping" procedure (adding or subtracting 
Y at discontinuity points) to enforce a smooth evolution oi (p{t). 

Using this result, the main axis vectors of the locally matching ellipse at 
instant t are obtained by going backwards along Xt{6) by an amount oi (p{t) or 
forwards by y — ^Pit)' i-e- 

a{t) =xt{-cp{t)) 

and 

If (p{t) has been adjusted for a smooth evolution over time, the same can be 
expected from the resulting a{t) and b{t). It is, however, not clear from this 
definition which one of these vectors specifies the major and minor axis of the 
ellipse, respectively, and it is possible that over the course of time the two vec- 
tors change roles. For a specific application of this result, further processing 
may therefore be necessary. 

The generalization of the instantaneous phase concept to the multivariate 
case may be complemented by the definition of a multivariate instantaneous 
amplitude A{t) such that 

z(f)=A(0 exp(i<^(0), 

which is given by 

A{t) =a{t)-ib{t). 

The channel-wise modulus of this quantity corresponds to the instantaneous 
amplitudes Aj(f ) of the component signals, while the argument comprises the 
phase differences between the global and the component signal instantaneous 

phases, (pi{t) — (p{t). 
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